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1.0 Problem Description; 

The Maypole Deployable Reflector is a large space structure 
consisting of a column and hoop structure which supports a circular 
antenna of 30-100 meters dia. It is a flexible structure and, hence, 
is described by partial, rather than ordinary, differential equations; 
such systems have a distributed parameter nature. 

Flexible structures like the above can be described by 
partial differential equations of the form: 


^t ^0\ 


+ AqU 


= F 


( 1 . 1 ) 


where u(x,t) is the displacement of the structure subject to external 
forces F(x,t). The structural stiffness is determined by the dif- 
ferential operator and associated boundary conditions. An example 
of such an operator for a membrane antenna is the following: 


AqU = V*E(x)Vu 


( 1 . 2 ) 


where E(x) is the distributed "stiffness" of the membrane and 
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Vu . gradient of u = . 


The damping in the structure occurs due to material properties and 







construction techniques. It is represented by the differential 
operator but the actual form of in a given application is 

much more difficult to determine than that of the stiffness A^. 
Some forms that might take are the following: 


DqU^ = (viscous damping) 

1/2 

< ^o'^t ~ ^0 ’^t damping) 


V. 


°0"t 




(?) 


Of course, there are other possible forms for and, in general, the 
damping operator might turn out not to be a differential operator at 
all; this is the case for some types of damping which are related to 
frequency in a very nonlinear way. 

The distributed parameter system (DPS) represented by (1.1) 
may be put into stat-e variable form: 


1^^ = Av(t) + Bf(t) 
. y(t) = Cv(t) 


(1.3) 


and B, C represent the control actuator and 
sensor influence operators, respectively, with f(t) being the vector 

of M actuator commands and y(t) the vector of P sensor outputs. It is 

assumed that the operator A in (1.3) generates a CQ-semigroup D(t) 

At 

which is similar to the matrix exponential e for lumped parameter 


where A = 


.-^0 -°o 


systems. 



2 . 0 Reduced-Order Modeling and Controller Design 


Many reduced-order modeling techniques exist for large- 
scale systems; of particular interest are the ones based on as}nnpto- 
tic methods, such as multiple time scales and singular perturba- 
tions. Let and be the subspaces of the total state space H 
with dim = N < « and H = ©■ Define the projection operators 

and (not necessarily orthogonal) and let 

V = P V. This decomposes v into v = + v_ and the system (1.3) 

K i\ N K 

into 
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y = s + ^R ^R 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


The reduced-order model (ROM) for the system is (2.1) and (2.2) with 
and assumed to be zero: 

r 

\ '*■ ®N ^ ^0 

(2.4) 

y 

Thus, the ROM depends on what choice of sub space is made and what . 
type of projection P^^ is used (or, alternatively, what is). The 
subspace H^j is called the ROM subspace and the subspace H^, the 
residuals subspace . The terms and A^^ v^ are called modeling 



error and B„ f and C„ v„ are called control and observation spillover, 
" ■ R R R 

respectively. This gives a general format for reduced-order modeling. 

The Controller Design is obtained by assuming the ROM (2.4) 
is the total system and using a linear controller of the form: 


f = Hii y -h ^ 
z = F z + Hy + Ef , z(0) = 0 


(2.5) 


where dim z = S N. 

All controller synthesis assumptions are made in terms of the 
ROM parameters (A^, Cj^, D) alone. This allows the synthesis to be 

carried out when these parameters are known. When the ROM parameters 
are not known or are poorly known, they must be estimated on-line or 
the controller must adapt itself in the presence of unknown parameters. 

3 . 0 Digital Parameter Estimation and Control 

The implementation of any control or parameter estimaticn 
scheme will be achieved with one or more on-line digital computers; 
consequently, all such schemes must be discrete- time. A discrete-time 
version of (1) may be obtained by assuming a constant input f (k) over 
the uniform time interval (k-1) At _< t < k At: 


v(k+l) =■ $ v(k) + Hf(k) 
y(k) == C v(k) 


(3.1) 


rAt 


where $ = U(At) and H = 


0 


U(t) Bdx. Other discrete-time versions 



can be obtained using non-uniform time steps. 


Now the model reduction described in (2.0) becomes: 

✓ 

v^Ck+l) = % Vj^(k) v^(k) + (k) 

< Vj^(k+1) = Vj^(k) + v^(k) + H^f(k) (3.2) 

y(k) = v^(k) + v^(k) 

\ 

where = P„ $ P„, etc. The reduced-order model is 
N N N 

v^(k+l) = v^(k) + H^f(k) 

(3.3) 

y(^y = S 

It is this approximation of the flexible structure whose parameters 
($^, Cj^) we will want to estimate on-line and use in any control 

schemes. The dimension N of this reduced-order model is related to 
the on-line computational capacity. 

Most parameter estimation schemes for lumped parameter sys- 
tems depend on an Autoregressive Moving Average (ARMA) model of the 
systeim. For distributed parameter systems, such as flexible struc- 
tures, we have developed the quasi-AKMA : 

N N 

y(k+N) = I a y(k+r-l) + J 3 f(k+r-l) + R(k) (3.4) 

r=l r=l 

where the coefficients a and 3 are obtained from (3.3) and the 

r r 



residual interaction term R(k) is given by 


N 

R(k) = C^Vj^(k+N) + I \ Vj^(k+r-l) 


( 3 . 5 ) 


with A = r - a C and 6 = T 

r r nr r r r r N 

By ignoring the term R(k) , we have an approximate ARMA model 
of the flexible structure whose parameters and 3^ we can estimate 
via least-squares or model reference techniques. However, the 
residual interaction terms will cause errors in these parameter esti- 
mates; hence, the effects of these terms must be analyzed to determine 
acceptable levels of residual interaction and tc. aid in the selection 
of the model dimension N. Compensation in the form of (digital) pre- 
filtering may be added to reduce the residual interactions with the 
parameter estimation scheme. 

Once the parameters are adequately estimated, they may be 
used in a discrete-time feedback controller of the form: 


f(k) = y(k) + 

z(k+l) = y(k) + H 22 z(k) 


( 3 . 6 ) 


Such a controller can substantially improve the performance of the 
flexible mechanical structure, as long as control or observation 
spillover is not excessive. Spillover is related to the model dimen- 
sion N and the actuator-sensor locations which are part of the control 


design. 



con- 


The use of adaptive structure controllers, i.e., 
trollers which tune themselves to the structure, raises even more 
complicated residual interaction issues. In particular, it has 
been shown via numerical simulation that adaptive cohtrol based on a 
reduced-order model can caiise unstable operation of the flexible 
structure if the spillover issue is not carefully considered. 

4 . 0 Recent Progress 

In the preliminary phase of this study, we have decided to 
consider two phases of operation of the Reflector: deployment and 

on-orbit operation. Deployment would be characterized by (possibly) 
high transient disturbances and rapidly-changing' configuration. 

Control applied during this phase would need to be done on a fast 
time-scale; hence, a low-order model of the reflector dynamics would 
be incorporated in the identification algorithm to enable the calcu- 
lations to proceed rapidly. On the other hand, once deployed the 
reflector would operate on-orbit in a relative steady-state. Hence, 
the identification (and control) could be done in a quasi-static mode. 

We feel this time-scale separation is a useful way to break 
down the identification problem. In the dynamic (deployment) mode, 
we have developed rapid, reduced-order identification schemes based 
on reduced-order lumped parameter models of the distributed parameter 
description of the reflector. I discussed some of the theoretical 
developments on this at the Yale Adaptive Control Workshop. 

For the time being, we consider the reflector to be a 
circular membrane; later, as the modeling aspect of the study develops. 



we expect to use more detailed models. Nevertheless, the generic 
behavior of a structure like the reflector can be revealed by even 
simple distributed models. In the future, I think we will try to 
fit distributed parameter models to the finite-element data on the 
reflector. 

4.1 Quasi-Static DPS 

In the quasi-static (on-orbit) mode, we do not have to rush 
to identify and control. Consequently, a least squares approach 
where a distributed parameter (partial differential equation) model 
is fitted to the sensor data seems reasonable. In discussions with 
H. T. Banks in May, we decided that a singular perturbations approach 
might yield some useful model reduction since the stiffness of the 
membrane is quite high along the radii where stringers are attached 
but much lower elsewhere on the reflector. We are using this idea to 
reduce the computational load of the identification problem. 

The quasi-static (or steady state) identification problem 
can be described by setting u(x,t) = u(x), i.e., no time-variation, 
in (1.1): 

Aq u = F = Bf (4.1) 

where Aq is described, for example, as the membrane operator (1.2) 
with E(x) unknown. Assume 



where are unknown scalars but are known functions (e.g., 

cubic splines or trigonometric functions) . This reduces the problem 
to the estimation of the parameters a 

From the Galerkin or finite-element method, we approximate 
u(x) in (4.1) by 


N 


(4.3) 


where 6|^(x) are known fimctions (although 0^^ is not necessarily equal 
to in (4.2)). The Galerkin method yields: for 2. = 1, 2, ... , N 


\ -=*0 \ 


where <u,v> = 


u(x) v(x) dx is the 




being the structure described in R" 
and (4.2) into (4.4) yields: 


2 

inner product on L (fi) with Q 
2 

(or R ). Substitution of (1.2) 
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a . 
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N 

k=l 




<F, e,> 


(4.5) 


where the known functions yield 
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j=l 


96j^(x) 3(}i^(x) 

9x . 9x , 

3 3 


+ <f>j^(x) A9j^(x) 


3 

Au = I 
j=l 



and 



From (4.3), it is possible to obtain the by knowledge 
of u(x) at various points on the structure (at least N of them) if 
we assume u(x) = u(x) ; this is true for N sufficiently large. There- 
fore solve: 

N 

I u^ ^^^i^ - u(x^) (4.6) 

k=l 

for u, when 1 < i < N. Furthermore, we need knowledge of the loads 
k — — 

f(x) applied by the actuators; this can be approximated and, hence, 
we assume <F,6j^> is also known. Then (4.5) takes the form: 

L a = F (4.7) 

where L is an N„ x N matrix and a, F are N„ x 1 and N x 1 vectors 

h — — £j 

respectively. Since in general N„ N there may be no unique solution 
for (4.7); however, a "least-squares" solution can be obtained via the 
pseudo-inverse matrix of L, i.e. , 

a = F (4.8) 

is the best mean-square solution of (4.7). 

The above describes our basic approach to the quasi-static 

A 

problem. Once the values of a, are determined from (4.8), they are 

used to approximate E(x) by a. ^.(x). Once E(x) Is identified, 

i=l ^ ^ 

it is possible to develop a quasi-static DPS control approach for (4.1). 


Thus, the control problem becomes the determination of actuator 



control commands f^, ... , (4.1) such that 

min II u(x) - u (x) ||^ 

f . 

1 

is achieved where u(x) is the actual antenna shape and is the 

desired shape, e.g., parabolic or spherical. We are in the process! 
of evaluating this identification and control approach for the 
circular membrane model of the antenna, llie singular perturbations 
idea for the reduction of (1.1) into two separate uniform problems 
is also under consideration. 

Results on the convergence and "well-posed"-ness of the 
above scheme and other similar ones are presented in "A Survey for 
Parameter Estimation and Optimal Control in Delay and Distributed 
Parameter Systems" by H. T. Banks, ICASE-NASA Langley Research Center, 
Report No. 81-26, August 17, 1981. 

4. 2 Dynamic DPS 

For the dynamic case (in which u(x,t) is a function of time, 
as well as x) , it is necessary to identify the terms in the linear 
differential equations that define the structural time variations in 
(1.1). In discussing potential identification procedures, one must 
keep in mind that identification of the dynamics of a large space 
structure is usually performed for one or more of the following reasons 

(1) To build a simulation model of the system which can 
be used for predicting the response to various 


types of inputs . 



(2) To build a controller for the structure. 

(3) To design adaptive control to account for structural parameter 
changes. (This purpose clearly requires an on-line 
identification procedure.) 

Relative to these purposes, various identification and 
adaptive control procedures have been considered. These include: 

(1) Adaptive observer - this procedure can be used either 
off-line or on-line. 

(2) Autoregressive moving -average (ARMA) identification - 
such an approach can be done off-line or recursively 

'I 

on-line. It has been applied by several investigators 
to adaptive control. 

(3) Frequency domain identification - although this is 
strictly an off-line procedure, it has the potential 
for producing a high order model without the inherent 
computational problems of many on-line procedures. 

(4) Indirect or implicit adaptive control - this procedure 
is useful for on-line direct computation of the con- 
trol gains without the need for explicit system para- 
meter identification. 

These procedures will now be discussed in more detail. 

(1) Autoregressive Moving Average (ARMA) Identification (Balas-Kaufman) 

A typical ARMA model between some output quantity y and a 
forcing function f can be written as: 

N N 

y(k+N) = I a. y(k+i-l) + I b. f(k+i-l) (4.9) 

i=l i=l ^ 


+ R(k) 




where R(k) Incorporates the residual effects. If R(k) can be modeled 
as the correlated noise sequence 

N 

R(k) = I c e(k+i-l) 
i=l 

where e^ is a sequence of independent Gaussian variables, then maximum 
likelihood estimation can be used. Thus, a., b., c. would be determined 

X I X 

so as to minimize 

L = 1/2 I (y(k) - i(k/k-l))^ 

A 

where y(k/k-l) is the predicted Kalman estimate of y. The optimum model 
order N can then be determined by evaluating Akaike’s criteria 
(- Log L + 2 * number of parameters) which weights both the error 
index and a measure of model complexity. 

C 2 ) Frequency D omain Pr ocedu re s (Kaufman) 

Alternately, if sinusoidal test signals can be used, it is 

suggested that frequency domain identification procedures also be 

1 2 

considered. Previous results ’ have shown that for flexible air- 
craft this approach is advantageous in that 

• the input design problem is eliminated; 

• fewer parameters need be identified per computational 
cycle; 

• Relatively simple algebraic least squares approaches 
can be used very effectively; 



• results from several controllers excited separately 
can be combined to improve accuracy. 

One disadvantage, however, is that the number of independent 
sensors must be equal to the assumed process order. Furthermore, this 
procedure is strictly an off-line type procedure and therefore is not 
suitable for use in an adaptive control scheme. 

( 3) Direct or Implicit Adaptive Control (Kaufman-Balas) 

An adaptive control procedure which appears appropriate for 

large space structures is the direct or implicit adaptive controller 

3 

proposed by Sobel, Kaufman, and Mabius. This, is a model reference 
procedure that adaptively tunes control gains so that the structural 
outputs follow the corresponding outputs of a reference model (which 
can be of lower order than the process). Stability is guaranteed 
provided that the structure has an input-output transfer function 
matrix that is positive real for some feedback gain matrix. Pre- 
liminary analysis indicates that this positive real property should hold 
for large space structural components provided that actuators and sen- 
sors are collocated and provided that velocity sensors are available. 

(4) Adaptive Observer (Balas) 

On-line adaptive parameter estimation of the DPS (1.3) cannot 
be done since this system is infinite dimensional. However, an adap- 
tive observer based on a reduced-order finite-dimensional model (2.4) 


can be accomplished. Such an adap t i ve ob se r ve r has the form: 





< 


+ bj^(t) f(t) + w(t) 


y(t) = c^ Vj^(t) 


(4.10) 


j ^ 

where K^ = [-k^ | Q^] is a known stable matrix and are the 

parameter estimates obtained from the following nonlinear adaptation 
laws ; 


a^(t) = (y(t)-y(t)) v^(t) 

bN(t) = (y(t)-y(t)) v^(t) 


(4.11) 


1 2 

where w(t) in. (4.10) and v (t), v (t) in (4.11) are auxiliary signals 
obtained by filtering the inputs f(t) and the outputs y(t). 

The parameter errors introduced by the reduced-order observer 
have been analyzed in Refs. 4 and 5. 

(5) Parameter Estimation Using a Linear Reinforcement Learning 
Factor (Desrochers) 

Recursive estimates which minimize 
f T ? / T " *v2 .N-1 

J = 2 (2^i ^i"^P “ 

(4.12) 

0 < a < 1 

S 


are considered where x^ and y“ are known and a^ are the parameters 



to be determined. This leads to the following recursive algorithm 


for parameter estimates : 

' ■ \ Vi ■ 

(xj V, - y*) (4.13) 

i 

P, = P, T P, T X, (1 + xJ P, T X, )"^ . 

k 1-a k-1 .2 k-1 k 1-a Tc k-1 k' 

(1-a) 

T 

^ \ \-l 

The effect of a (Forgetting Factor) is to weight newer data more 
heavily than past measurements in order to better adapt to changing 
parmaeters. 

Of course, this approach has the disadvantage that the 
effect of noise will also be detected and used to modify the estimate 
a. Therefore, the scheme functions best when the parameter variations 
are larger than the residual fluctuations due to noise. The factor a 
is a variable and it is updated through a linear reinforcement algo- 
rithm . Such an algorithm makes decisions about the reliability of the 
present data and adjusts the factor a accordingly. 

(6) Feedback Control (Balas) 

Analysis of the stability of feedback control based on 
finite element (Galerkin) reduced-order models of DPS such as flexible 
structures is considered in Ref. 6. This analysis assumes that the 
identification problem has been successfully completed by one of the 
above-mentioned procedures. Control procedures for DPS are evaluated 


in Ref. 7. 




(7) Reduced-Order Modeling for Nonlinear Systems (Desrochers) 


The first case is concerned with nonlinear systems that 
can be modeled by 

n 

x(k+l) = I A. F,(x.(k)) + Bu(k) (4.14) 

j^l J 2 2 

where Aj is an n x m matrix of constants, x(k)eR^, (•) are real 
valued nonlinear vector functions of dimension m, B is n x 1, and u(k) 
is the scalar input. Consider the equation error for a single non- 
linearity in (4,14) 


e.(k) = x(k+l) - [A. B] 
2 2 


F^(k) 

u(l) 


j = 1, 2, 


n 


(4.15) 


Minimizing 


N-1 

V 

L 

k=0 


eT(k)e^ (k) 


(4.16) 


leads to an algorithm for assigning costs to the Fj(’) order to 

retain the most dominant non-linearities. Since each F.(') is a 

3 

function of only one x. this leads to the retention of the most domi- 

3 

nant states. 

The more general case of 


n 

x(k+l) = ^ A F (x ,X2> ... , X ) + Bu(k) 

3=1 ^ ^ 


(4.17) 



can also be handled this way. Note that now, model reduction may 
not necessarily lead to an elimination of state variables, but the 
number of terms, n, will be reduced. 

These techniques are described in Refs. 8, 9 and 10. We 
feel that in the case of nonlinear deformation of the antenna it 
is wise to make a nonlinear model reduction before linearizing the 
system for control. Also, the nonlinear reduction can be used to 
obtain a nonlinear version of the ARMA described above, this could 
be useful in the identification of nonlinear models of the DPS. 
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